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Aperiodic dynamical decoupling (DD) sequences of vr pulses are of great interest to decoherence 
control and have been recently extended from singie-qubit to two-qubit systems. If the environmen- 
tal noise power spectrum is made available, then one may further optimize aperiodic DD sequences 
to reach higher efficiency of decoherence suppression than known universal schemes. This possi- 
bility is investigated in this work for the protection of two-qubit states, using an exactly solvable 
pure dephasing model including both local and nonlocal noise. The performance of optimized DD 
sequences in protecting two-qubit states is compared with that achieved by nested Uhrig's DD 
(nested-UDD) sequences, for several different types of noise spectrum. Except for cases with noise 
spectrum decaying slowly in the high-frequency regime, optimized DD sequences with tens of control 
pulses can perform orders of magnitude better than that of nested-UDD. A two-qubit system with 
highly unbalanced local noise is also examined to shed more light on a recent experiment. Possible 
experiments that may be motivated by this work are discussed. 

PACS numbers: 03.67.Pp, 03.65.Yz, 07.05.Dz, 02.70.-c 



I. INTRODUCTION 



In one way or another any two-level system (qubit) 
is coupled to some environmental degrees of freedom. 
This inevitable system-environment coupling leads to de- 
coherence. To protect quantum coherence as a key re- 
source for quantum technologies, many schemes have 
been proposed to dynamically eliminate the unwanted 
qubit-environment coupling [l|-l5|. 

Analogous to the spin-echo technique widely adopted 
in nuclear-magnetic-resonance studies @, various dy- 
namical decoupliirg (DD) sequences of instantaneous con- 
trol pulses 0, 0-[l| have been extensively studied. In 
particular, due to Uhrig's DD (UDD) sequence ^, B , re- 
search activities focusing on DD have surged recently [lOl - 
[l5j . In addition to its high efficiency in theory, the power 
of UDD lies in its universality [l6l . \n\ . Indeed, the work- 
ing mechanism of UDD does not rely on detailed assump- 
tions about the system-environment coupling or about 
the environment. Nevertheless, if the actual form of the 
noise spectrum of the environment is available, then a 
locally optimized DD (LODD) sequence based on 
the noise spectrum can outperform UDD, with the pulse 
locations optimized according to an exact decoherence 
function |18l - [2(]| . The reason for the success of optimiza- 
tion is simple. In suppressing the pure-dephasing of a 
qubit, a UDD sequence minimizes a decoherence filter 
function (defined later) in the neighborhood of zero fre- 
quency but gradually becomes less effective for large fre- 



quencies. As such, if the noise spectrum of the environ- 
ment is known, then its actual behavior at appreciably 
nonzero frequencies makes room for further optimization 
of DD sequences. This optimization approach is some- 
what in the same spirit of the continuous DD approach 
d, 0] , insofar as both attempt to make the full use of the 
noise spectrum. 

Extension of DD to two-qubit (or even multi-qubit) 
systems is crucial towards efficient protection of quantum 
entanglement. For a known initial state of a two-qubit 
system, an extended UDD sequence is found by involv- 
ing nonlocal control operators, with its theoretical perfor- 
mance essentially identical with that in one-qubit UDD 
For general situations, nested-UDD sequences. 
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cases 

initially proposed for suppressing both dephasing and re- 
laxation in one-qubit systems [22| , are advocated to pro- 
tect two-qubit states in a universal manner, i.e., without 
knowing the details of the system-environment coupling 
or of the noise spectrum [2j, [23|. For example, it was 
shown that to lock an unknown superposition state of 
two known basis states to the A^th order, three layers 
of UDD sequences and hence about control pulses in 
total will be needed [2l[. If the initial state is totally 
unknown, then to reach the same level of decoherence 
suppression one needs four layers of UDD sequences and 
hence about N"^ control pulses [1^ [l^l • One important 
question then arises. That is, if the noise spectrum of 
the environment of a two-qubit system is available, then 
can we significantly improve entanglement protection by 
further optimizing the locations of the instantaneous con- 
trol pulses as what was done in [3 single-qubit 
systems? If yes, then the required number of pulses can 
be much less and more understandings of entanglement 
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protection might emerge. Using an approach extended 
from the above-mentioned LODD for single-qubit sys- 
tems, this question is answered here via a pure dephasing 
model of an open two-qubit system. The performance of 
optimized DD sequences in preserving two-qubit states 
is compared with that achieved by nested-UDD, for sev- 
eral different types of noise spectrum. Except for the 
Lorentzian type of noise spectrum, it is found that opti- 
mized DD sequences can protect two-qubit states orders 
of magnitude better than nested-UDD. As a result, on 
the one hand nested-UDD can be seen as a powerful and 
a general-purpose scheme for protecting two-qubit states, 
and on the other hand the optimized DD approach can 
be seen as system-specific DD schemes with even better 
performance. In addition, to shed more light on a re- 
cent experiment of entanglement protection via DD [2^ , 
a two-qubit system that shares important noise features 
with the experiment is studied. 

Our plan for this paper is as follows. In the next sec- 
tion, we discuss a model that describes pure dephasing 
processes of a two-qubit system in the presence of instan- 
taneous 77 pulses. Based on exact expressions of decoher- 
ence effects, we elaborate in Sec. Ill our optimization 
procedure to find the optimized pulse locations. In Sec. 
IV, the performance of optimized two-qubit DD sequence 
is illustrated in four subsections treating different types 
of noise spectrum. Sec. V concludes this work and pro- 
poses two types of experiments. 

II. PURE DEPHASING MODEL OF 
TWO-QUBIT SYSTEMS 

In general, decoherence as a complicated process in- 
volves both population relaxation and dephasing. Yet, 
if we increase the strength of an external polarization 
field such that all the energy level splittings are suffi- 
ciently large, then the relaxation can be made negligible 
within a time scale of interest and dephasing becomes the 
only source of decoherence. Under such a pure dephas- 
ing assumption, the environmental noise may be modeled 
as classical random fields causing random phase shifts, 
a valid treatment in many depha sing environments like 
spin bath in solid-state systems [13, [13, Udl ^^'^ back- 
ground noise in superconducting qubits [26[. Here we 
follow the methodology proposed in [2^ . We can then 
write the pure-dephasing Hamiltonian of an open two- 
qubit system as 

H - fi{t)a,, + f2{t)(T,, + h{t)<j,,a,,, (1) 

where fi (t) are assumed to be independent random noise 
variables with a Gaussian distribution. We further as- 
sume 

(/.) - 0, (2) 

{f^itl)Mt2))=g^itl-t2), (3) 

with the correlation function gi{t) being an even func- 
tion. The fi{t) [,f2{t)] term in the Hamiltonian depicts 



the noise seen by the first (second) qubit alone; whereas 
the fslt) term reflects how noise may jointly impact the 
two qubits. In the following we loosely call the f3{t) noise 
term as a term of nonlocal noise. If the two qubits are 
far apart, then this nonlocal noise term should be very 
small and can be neglected. As a consequence the deco- 
herence control problem is expected to share important 
features with single-qubit DD (note that fi{t) and 
f2{t) are assumed to be independent here). However, of 
our concern in this study are two nearby qubits inter- 
acting with each other, and therefore the nonlocal noise 
term should be included to account for random fluctua- 
tions in the qubit-qubit mutual interaction (for example, 
due to the fluctuations in qubit-qubit distance caused by 
lattice vibrations, or in the context of super-conducting 
qubits ^2^ , due to the fluctuations in a third device that 
is responsible for qubit-qubit coupling). Certainly, in a 
more realistic environment the fsit) term might be cor- 
related with the local noise terms. Such kind of potential 
correlations are neglected in our model. From a differ- 
ent perspective, one may regard the first qubit as a part 
of the environment of the second qubit, and then the 
/a {t) term models how the dynamics of one qubit might 
change the environment experienced by the other qubit. 
Throughout we assume dimensionless units. 

Let I t) and | |) be the eigenstates of az , then a general 
two-qubit pure state at time zero can be written as 

|f (0)) = a\ U)+/3\ ;t) +7l n)+v\ n), (4) 

where the upward or downward arrows in each of the four 
components represent the spin states of the two qubits. 
The Hamiltonian in Eq. (1) then gives rise to the follow- 
ing state vector at time t, 

+ ,,e-i[^^(*)+^^(*)+^^«l|n), (5) 

where Fi{t) = fi{t )dt . To analyze the density matrix 
p{t) = |4'(i)}(^(t)| averaged over noise histories (denoted 
(•)), we further define four basis states |0} = | |i), |1) = 
It) J2) = I ti), |3) = I tt)- Then ah the averaged den- 
sity matrix elements (in the absence of DD control pulses) 
can be easily worked out. For example, the mean value 
of pm{t), denoted poi(0, '^a-'^ expressed as 

poi(t) = «*^(e-i[^^^(*)-^'3(t)]g-i[J^2(0-F3(t)]^ 
= a*/3e-2([^'2(t)-i=^3(t)]^) 
= a*/3e-2[<^l(*))+<^l(*))l. (6) 

Here, in obtaining the last equality we have used the 
relation (/2/3) = as well as the Gaussian nature of the 
noise. Similar expressions can be obtained for all other 
density matrix elements. 
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In our pure-dephasing model, the error operator for the 
first (second) qubit is a^^ (czs) only, whose detrimental 
effect can be suppressed by a local control operator 
{(7x2) ■ As seen from Refs. |23l. [23| ■ a nested- UDD scheme 
with two layers of CTxi and pulses can suppress the de- 
phasing to the A^th order, with about iV^ pulses in total. 
To examine if we can further improve the performance 
by optimizing the pulse locations, let us now consider 
the following scenario: n pulses of tt rotation along axis- 
X are applied to the first qubit, with the pulse locations 
given hy ti,t2, ■■■,tn, whereas m analogous pulses are ap- 
plied to the second qubit at times tj^' , , ^m' • As a 
result, totally n + m pulses are applied to the two-qubit 
system. We arrange the n + m pulse timings in increas- 
ing order and denote them by t-^" , 



t" < t 



(i+i)" • 



.., t, , ■," , with 
At each of such instants, either the 



or operator switches its sign. At the same time, the 
operator (Tz-^crz2 changes its sign n + m times at these 



instants, 
tions [^]: 



This motivates us to define three switch func- 



si(/) tfci < i'<i(fei+i), fci = 0,l,...,n, 

S2(i') = <i'<t(fe,+i)', fc2-0,l,...,m, 

S3(i') = ifcf; <t'<t(fc3+i)", fca = 0,l,...,n + m, 

with to = 0, and t^+i = = t{n+m+i)" ^oi 

times outside the domain [0,t\ these switch functions are 
defined to be zero. The influence of the DD pulses can 
then be expressed in a rather compact form. Still taking 
the decay of poi as an example, we obtain 



PoiW=a*/3e^2[(F|(t))+(^|(t))]^ 
where Fi{t) is given by 

F^{t) = / Mt)s,{t)dt. 



(7) 



(8) 



To proceed we next deflne three fllter functions [9| from 
the Fourier transform of Si{t), i.e.. 



(9) 
(10) 



S3(i')e'"''di' = -yn+rniujt). (11) 



Here, the filter function with M (which is n, m, or n + m) 
pulses is defined as 

M 

VMii^t) = 1 + (-l)^'^+iei"* + 2^(-l)^ei"*^% (12) 

with Sj being the pulse location normalized to the to- 
tal duration t, i.e., Sj = tj/t. The decay factor in the 
exponential of Eq. ([7]) then reduces to 

{Fiit)) + {Fiit)} 

dtidt2 32(^1)32(^1 -t2)s2{t2) 



00 pOO 



00 J —00 
00 pOO 



00 J —00 



dtxdt^ 53(^1)53(^1 -^2)53(^2) 



\y,j,(uji)\ — 



1 



^1 I / \ |2 '^3(^) 

-f- / |?/,„+„(cji)| 2 







(13) 



where Si{iS) is the noise spectrum (or spectral density), 
namely, the Fourier transform of the noise correlation 
function 5i(t), with 



9i{t) ^ — I Si{uj) cos{ujt)duj. 



(14) 



Note that an equivalent expression may be derived us- 
ing an exactly solvable spin-boson model with pure- 
dephasing. Interestingly, the decay factor in Eq. (O is 
seen to be the sum of two decoherence functions, each 
being analogous to that in the previous DD studies ^ ^] . 
Yet, the problem here is still rather complicated for two 
reasons. First, for the same control pulse, the coefficient 
(—1)'' can take different values in the three filter func- 
tions. Indeed, its actual value depends on its relative 
pulse location in the time sequences of tk^ , of i^' , or of 
t J," . Second, density matrix elements do not decay in 
the same fashion, and hence the decay of all off-diagonal 
density matrix elements needs to be accounted for. 

Repeating the above procedure for all other density 
matrix elements, we finally obtain the averaged full den- 
sity matrix p{t) as follows: 



J 



/ |a|2 a*/3e-r2(*)-r3(t) 

/3*ae-r2(t)-r3(t) \p\2 

^ 77*ae-r2(t)-ri(t) ,^*^g-ri(t)-r3(t) 



a*"fe 


~ri(t)-r3(t) 


a*r]e~ 


-r2(t) 


-ri(t 




-r2(t)-ri(t) 




-ri(t) 


-r3(t 








-r2(t) 


-r3(i 




-T2{t)-T^(t) 
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with the decay exponents for the off-diagonal elements 
given by 



r2 = 



\yn\ —a^. 



UJ'- 
I2 5'2(w) 



(kjj. 



(15) 
(16) 
(17) 



For convenience, unimportant constants in Eqs. (|15p - p7)) 
are either set to unity (e.g., t = 1) or absorbed into the 
noise spectrum. 

Before ending this section, we stress that our straight- 
forward calculations above are made possible by first as- 
suming the statistical independence of /i {t) and /2 {t) . If 
/i (t) and /2 (t) has nonzero correlations (noise under this 
correlated situation may be also called nonlocal noise, 
which is much different from our case here), then many 
of the density matrix elements cannot be evaluated ana- 
lytically. 



III. OPTIMIZATION PROCEDURE 

Next we aim to optimize the pulse locations to keep 
p{t) as close as possible to the initial state. Taking 
the trace fidelity C{t) — Tr[p(t)p(0)] as a measure of 
DD performance, it is straightforward to carry out the 
optimization if the initial state is known. However, in 
many cases a two-qubit state to be protected or stored 
is unknown, and as such the average fidelity for all pos- 
sible initial states can be of more interest. Averaging 
C{t) over all initially pure states and using the fact that 



we arrive at 



(18) 



The optimization then becomes the minimization of 
$(t) = 4[1 - C{t)\. The function $(i), called the per- 
formance function below, is given by 



$(i) = 3 - (e" 



-r2 



(19) 



Clearly < ^{t) < 3, and a smaller value of ^{t) indi- 
cates a higher degree of decoherence suppression. There- 
fore $(<:) can be also regarded an error function. 

To minimize the value of the performance function, we 
first consider a total of n -I- m pulses, among which m 
pulses are applied to the second qubit. Assuming that 
the n + m pulses are applied at the timings 



Si < 62 < Ss < ... < 6, 



one first needs to pick out m pulse locations to apply 
the control operator. The number of such choices is 
given by the combination C^„. Then the analytic form 
of the three filter functions ?/„(a;), ym(j-o), and y„+m('^) 



TABLE I: Impact of nonlocal noise on the performance of 
optimized two-qubit DD sequences. The local noise spectrum 
is assume to be Si{lu) = ajO(l — cj) for the first qubit and 
S2{ijj) = U)Q{1 — oi) for the second qubit. n + m = 8. In 
the column of pulse location, a listed integer j means that 
the jih pulse will be a pulse (i.e., applied to the second 
qubit). Performance is optimized by considering all C| and 
Cg possible pulse allocations for the two qubits. 



Nonlocal spectrum 53(1^) 
2aje(2 - uj) 
0.5cje(0.5 - oj) 
0.1a;e(0.1 - a;) 
-(infinite cutoff) 



Pulse location Performance 



no nonlocal noise 



2,4,6,8 
2,4,6,8 
2,4,6,8 
2,4,6,8 
2,4,6,8 



4.59 X 10"'^ 
4.59 X 10~^ 
4.43 X 10"^ 
1.67 X 10"^ 
4.08 X 10~^° 



can be fixed. We next minimize the performance func- 
tion $(t) by optimizing the (n + m) pulse locations using 
the line-search algorithm. In every iteration, the algo- 
rithm searches for the minimum along the direction set 
by a gradient. We use n + m equally spacing pulses and 
a two-layer nested-UDD sequence as our initial guesses. 
Once the result converges, we obtain the locally opti- 
mized locations for n + m control pulses in a particular 
(n -|- TO, to) scheme. To find the optimal pulse locations, 
we test all C^+m possible choices of pulse allocations for 
the second qubit. Furthermore, by scanning the value of 
TO from 1 to the total number of pulses, we can find an 
optimal pulse-number partition for a fixed total number 
of pulses. 



IV. PERFORMANCE OF OPTIMIZED 
SEQUENCES 

In the following, for convenience we set t — 1 when 
comparing cases of different noise spectrum. Effects of 
varying t can be understood as the result of a rescaling 
of the strength and shape of the noise spectrum [23| • To 
appreciate that the issue of two-qubit decoherence sup- 
pression is significantly different from a one-qubit prob- 
lem, let us first illustrate the influence of the nonlocal 
noise term /3(i)o'2iCrz2 on the performance of DD with 
n + m — 8. Computational examples are shown in Ta- 
ble HI It is seen that in the absence of nonlocal noise, the 
optimized performance (error) is of the order of 10"^", 
and the optimized DD pulses are applied at 

[0.10 0.10 0.35 0.35 0.65 0.65 0.90 0.90], 

which is simply two optimized 4-pulse sequences simul- 
taneously applied to the two qubits. However, this DD 
sequence cannot suppress any nonlocal noise because ev- 
ery pair of such simultaneous pulses will keep the sign of 
o'ziCr22 unchanged, and the associated decay exponent T3 
would be the same as that without DD, which is given 
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TABLE II: Performance of optimized two-qubit DD sequences 
as compared with nested-UDD, for a few examples of local 
and nonlocal noise spectrum. Performance is optimized for 
each combination taking into account all possibilities 

of allocating m pulses to the second qubit. In the column of 
pulse location, a listed integer j means that the jth pulse is 
a pulse (i.e., applied the second qubit) during the first 
half of the dynamics. The second half of the pulse sequence is 
symmetric to the first half. Ohmic noise spectrum with hard 
cutoff is considered. 



Si = ujO{1 — uj] 


),S2=ujQ{l~uj),S'i 


= 2tje(2 - cj) 


Pulse combination 


Pulse location 


Performance 


nested-UDD (2) 


3 


7.32 X 10"" 


Ci 


3 


8.66 X 10"^ 


cl 


2,4 


4.59 X 10"^ 


nested-UDD (3) 


4,8 


2.45 X 10"*' 


ch 


4,8 


3.04 X 10"'^ 


Cl, 


2,5,8 


6.14 X 10"^ 


Cfs 


1,3,5,7,8 


1.17 X 10"" 


Si = aje(l -a;),S2 = Lue{l - w), S3 = 


0.5a;e(0.5 - a;) 


Pulse combination 


Pulse location 


Performance 


nested-UDD (2) 


3 


3.26 X 10-* 


Ci 


3 


8.14 X 10"^ 


Ci 


2,4 


4.59 X 10"^ 


nested-UDD (3) 


4,8 


1.66 X lO"'' 




3,8 


1.88 X 10"'^ 


ch 


3,5,8 


7.06 X 10"" 


Cl, 


1,3,4,6,8 


6.26 X 10""^ 



by 

r,.4r%M5Hl>/a.„^ (20, 

Interestingly, as shown in Table HI even with a very weak 
nonlocal noise added, e.g., 53(0;) = O.lcij0(O.l— w), where 
O(-) is the Heaviside step function, the minimal error 
that can be achieved by an optimized DD sequence is al- 
ready increased by five orders of magnitude! This clearly 
addresses the importance of taking nonlocal noise into 
consideration for decoherence suppression. With this un- 
derstanding we are now ready to examine the optimiza- 
tion of DD in two-qubit systems. 

A. Spectrum with Hard Cutoff 

Previous work showed that in single-qubit cases, an 
optimized DD sequence can greatly outperform UDD for 
an Ohmic noise spectrum with hard cutoff [l^, [l8l - l20| . 
It is observed that this finding also holds in our two- 
qubit dephasing model. Here, we set the local noise 
spectrum of the two qubits to be the same, and then 
vary the intensity of the nonlocal noise spectrum, from 
twice as much as the local noise spectrum to a relatively 



TABLE III: Same as in Table II, but for Ohmic noise spectrum 
with hard cutoff at larger frequencies and for 1// spectrum 
with hard cutoff. 



Si = a;0(5 — Lo 


), 5*2 = a;B(5 — uj] 


],S3 = 


= a;e(3-a;) 


Pulse combination 


Pulse location 




Performance 


nested-UDD (2) 


3 




1.55 


Cl 


3 




0.80 


Ct 


2,4 




0.54 


nested-UDD (3) 


4,8 




0.36 




3,8 




6.63 X 10"^ 


ch 


2,4,6,8 




1.48 X 10"^ 


Si = ie(10-a; 


),52 = ie(10-a;),S3 


= ie(5-a;) 


Pulse combination 


Pulse location 




Performance 


nested-UDD (2) 


3 




0.61 


Cl 


3 




0.60 


Ct 


2,4 




0.41 


nested-UDD (3) 


4,8 




0.32 


ch 


4,8 




0.22 


ch 


2,4,6,8 




9.96 X 10"^ 



weak one. As seen from Table |TT1 the optimized DD se- 
quence is in general much better than nested-UDD. In 
particular, for a 15-pulse sequence (n + m = 15), the op- 
timization improves the performance by many orders of 
magnitude as compared with a two-layer nested-UDD(3) 
(here nested-UDD(fc) means that k pulses applied to the 
second qubit in the outside layer, and within each layer a 
/j-pulse UDD sequence is applied to the first qubit). As a 
matter of fact, by comparing Table HIl and Table IVll it is 
seen that the best performance of optimized 15-pulse DD 
sequence for 6*3 — O.5aj0(O.5 — uj) can even be two orders 
of magnitude better than nested- UDD (4), which requires 
24 pulses. This indicates that if the noise spectrum is 
known, then the required pulse number to achieve a given 
degree of decoherence suppression does not have to scale 
with N"^ as in nested-UDD(Af). Further, as suggested by 
the results in Table |lTl to achieve the best performance 
the number of pulses is larger than that suggested by 
nested-UDD. But interestingly, this does not necessarily 
mean that the number of pulses should be as close as 
possible to the number of Gx^ pulses. For example, the 
C^5 (instead of Cjg) case shown in the bottom of Table II 
gives the best performance. This is one of the subtle con- 
sequences of nonlocal noise spectrum. Lastly, as 6*3(0;) 
changes from 2a;8(2 — uj) to 0.5wO(0.5 — uj), it is seen 
from Table II that not only the optimized performance 
changes, but also the whereabouts of ax2 pulses become 
much different. The physical intuition why the cTx2 pulses 
should be shifted in this manner is far from obvious, re- 
flecting the importance of the relative locations of ax2 
pulses with respect to Gx^ pulses. 

In single-qubit cases [9l [l8l - [20| . LODD is more power- 
ful for a noise spectrum with a larger cutoff. This mo- 
tivated us to investigate some representatives cases with 
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0.2 0.4 0.6 0.£ 
nested-UDD(2) 



0.2 0.4 0.6 0.£ 
nested-UDD{3) 



0.2 0.4 0.6 O.S 
optimized (8,2) 



0.2 0.4 0.6 O.S 
optimized (15,3) 



0.2 0.4 0.6 
optimized (8,4) 



0.2 0.4 0.6 0.8 
optimized (15,7) 



FIG. 1: Comparsion of specific pulse locations between 
nested-UDD and our optimized two-qubit DD sequences. 
Si = uje{5 - uj), S2 = uje{5 - uj) and 5*3 = a;e(3 - cj). 
The left panel is for n + m — 8, and the right panel is for 
n + m = 15. The grey (red) lines stand for timings of the 
pulses applied to the second qubit. 



a larger frequency cutoff. Here we consider Ohmic and 
1/f spectrum in Table Hill fwith cutoff). Here the 1//- 
type noise is of interest because it has been observed in 
many solid state implementations, especially in super- 
conducting qubits [2a,[H,[33|- As seen from Table iHll in 
both cases the two-qubit dephasing can be greatly sup- 
pressed. Interestingly, although the Ohmic and the 1/f 
cases represent drastically different noise spectrum, in 
these two cases the optimized 15-pulse two-qubit DD se- 
quences give similar performance that is five or six orders 
of magnitude better than nested-UDD (3). Clearly then, 
the improvement afforded by optimization is sensitive to 
the cutoff frequency value, but rather insensitive to the 
shape of the noise spectrum before the cutoff. Similar im- 
provements have been observed previously in single qubit 
LODD 



In Fig. [l]we show the optimized pulse locations in the 
case of Si = uj&{5 — uj), S2 — w6(5 — w) and ^3 = 
w8(3 — u). In the case of n + m = 8 with rn = 2, the 
locations of all pulses in the optimized DD sequence 
differ remarkably from that in nested-UDD(2). In the 
case of n -t- m = 15 with m = 3, the first pulse is 
the third pulse in the optimized two-qubit sequence but 
the fourth pulse in nested-UDD(3). These observations 
might be relevant to future experiments. 



B. Spectrum with Soft Cutoff 

Noise spectrum with soft cutoff only gradually decays 
to zero. The associated task of decoherence suppression is 
more challenging. References [H, [H, showed that 
for single-qubit dephasing caused by soft-cutoff noise, 
UDD can only give a performance analogous to the Carr- 
Purcell-Meiboom-Gill (CPMG) sequence [3l|. Mathe- 
matically, this can be understood from the decay ex- 
ponents Ti determined by an integral of a filter func- 
tion multiplied by a noise spectrum. Because the noise 
spectrum is not rapidly vanishing in the high frequency 
regime, the behavior of a filter function in the high fre- 
quency regime may be important and as a result the op- 
timization becomes less effective. 

Here we study the optimization of DD using vari- 
ous noise spectrum combinations, involving both super- 
Ohmic (aj",a > 1) spectrum and Lorentzian spec- 
trum with soft cutoff. The results are presented in Ta- 
ble llVl In the first super-Ohmic case with an exponential 
soft cutoff, the improvement of optimized DD sequence 
over nested-UDD is still magnificent (note that a non- 
optimized DD sequence might even speed up the deco- 
herence process [S^l). Indeed, the exponential cutoff can 
be regarded as an intermediate case between a truely slow 
cutoff (e.g., Lorentzian spectrum) and a hard cutoff stud- 
ied above. The second case [Si{uj) = S2{uj) — ujQ{1~uj), 
S3 = Ji'^i ] should be compared with Table HIl with the 
only difference being that a nonlocal Ohmic spectrum 
replaced by a weak nonlocal Lorentzian spectrum. Con- 
sequences of this change in the nonlocal noise spectrum 
are: (i) Performance of nested-UDD is greatly reduced, 
so is the performance of optimized DD, (ii) With the 
same value of n-|-m, optimized DD can only improve the 
performance over nested-UDD slightly, (iii) Finally, the 
performance still improves with the increase of n+m, but 
very slowly. Similar observations can be made from the 
third case in Table IIV[ where the local noise spectrum is 
assumed to be Lorentzian. 



C. Asymmetric Local Spectrum 

So far local noise spectrum is assumed to be the same 
for the two qubits. In practice the local noise spectrum of 
individual qubits can be much different. The recent ex- 
periment on DD of a two-qubit system represents one ex- 
cellent example of unbalanced local noise spectrum [25i |. 
In the experiment the dephasing time of a nuclear spin is 
two orders of magnitude larger than that of an electron 
spin. The DD pulses are then applied to the electron 
spin only. After a two-flip control sequence, the lifetime 
of a pseudo-entangled state in the experiment becomes 
analogous to that of the nuclear spin, which implies that 
control of the nuclear spin is needed to further improve 
the entanglement protection. 

To mimic an environment like the experimental setup 
in Ref. ^5|], we consider an environment with strongly 
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TABLE IV: Same as 


in Table II, but involving 


; noise spectrum 


without hard cutoff. 








2 o 2 2 

,52 = a;3e-" , 53 = a;e-'" 


Pulse combination 


Pulse location 


Performance 


nested-UDD(2) 


3 


5.31 X 10"^ 


Ci 


2,4 


1.04 X 10"^ 


nested-UDD(3) 


4,8 


1.44 X 10"" 


c!s 


2,5,8 


5.25 X lO"'' 


Si = ioQil - 


~ tj), S2 = a;9(l - a;), 53 


0.2 


Pulse combination 


Pulse location 


Performance 


nested-UDD(2) 


3 


4.36 X 10"^ 


Cl 


2,4 


1.67 X 10"^ 


nested-UDD(3) 


4,8 


1.20 X 10"^ 


•-^15 


2,4,5,7,8 


4.74 X IQ-'* 


o 0.2 






Pulse combination 


Pulse location 


Performance 


nested-UDD(2) 


3 


2.87 X 10"^ 


Cl 


2,4 


2.08 X 10"^ 


nested-UDD(3) 


4,8 


1.36 X 10"^ 


ch 


2,4,6,8 


3.96 X 10"^ 



asymmetric local noise. For convenience the noise spec- 
trum is assumed to be a constant subject to a hard cutoff, 
i.e., Si{uj) = 109(10 - w), S'2(w) = 0.19(0.1 - uj), and 
53 (w) — 0.059(0.05 — ri). For these noise parameters the 
dimensionless dephasing rate of the first qubit is about 
two orders of magnitude larger than that of the second 
qubit. The optimization results are shown in Table IVl At 
least two interesting observations can be made. First, in 
the case of n+m = 4, assigning some control pulses to the 
second qubit may bring down the performance. Hence, 
if the local noise spectrum is strongly unbalanced and if 
n+m is highly limited, then it is beneficial to control only 
one of the two qubits. This is consistent with the exper- 
imental implementation in Ref. [25J. Second, when the 
minimized total error approaches the error caused by the 
free decay of the second qubit, solely applying pulses 
can no longer improve the DD performance. As shown in 
Table |Vl the DD performance without any pulses is 
bounded by about 10~^. However, if now two pulses 
are applied, then the DD performance can be improved 
by many orders of magnitude (e.g., in the 12-pulse case). 
Note also that for fixed n + m, increasing m too much 
undermines the performance. Again, this is because the 
number of pulses applied to the two qubits should be 
balanced so that errors of the two qubits are suppressed 
to the same level. Less pulses should be applied to the 
second qubit because it is more weakly coupled to its en- 
vironment. As such, the number of ax-^ pulses (n) should 
always dominate. In Fig. [5] we display the optimized 
pulse locations in several cases. The left panels of Fig. [5] 
represent optimized pulse locations if m — 0. The right 
panels show the optimized pulse locations if m = 2. In- 



TABLE V; Same as in Table II, but for an example where the 
local environments of the two qubits are highly unbalanced. 

Si = 109(10 - g;), ^2 = 0.19(0.1 - a;), ^3 = 0.059(0.05 - a;) 
Pulse combination Pulse location Performance 



Cl 


no pulse 


1.30 


ci 


2 


2.00 


ci 


no pulse 


2.00 X 10"^ 


Ci 


3 


7.64 X 10"^ 


ci 


1,3 


1.30 


ci 


1,2,4 


2.00 


C12 


no pulse 


1.99 X 10"^ 


C12 


4 


1.57 X 10"^ 


C12 


3,5 


6.25 X lO"'' 


<-^12 


2,4,6 


7.45 X 10"^ 


^8 

'-^12 


1,3,4,6 


1.29 



terestingly, in the upper right panel, it is found that for 
n + m = 4, the two pulses for the second qubit should 
be applied at the same time, thus canceling themselves. 
Our optimization procedure therefore requests a null ac- 
tion on the second qubit if n -I- m is small. For the case of 
optimized (8,2) shown in Fig. 2, the pulse combination is 
the same as in nested-UDD(2), but with optimized pulse 
timings. We have also compared the optimized (12,2) 
case in Fig. 2 with an asymmetric nested-UDD of total 
11 pulses, namely, UDD(3) on the first qubit in the inner 
layer and UDD(2) on the second qubit in the outer layer. 
The performance of this asymmetric nested-UDD is only 
0.517, still far from what is achieved here (only at the 
cost of one additional pulse in the middle interval of the 
inner layer). 



D. Limitations on Optimized DD 

As in single-qubit DD, we would naturally attempt to 
increase the total pulse number n -f m to achieve better 
and better DD performance. However, the optimization 
procedure faces much difficulty when n + m increases 
considerably. First, the optimization algorithm becomes 
much slower in each run and may fail to give a con- 
verged result. For example, for n + m — 24, the con- 
verged pulse locations might be inappropriate, i.e., vio- 
lating the initial ordering of the pulse locations. Worse 
still, for n + m = 24, we have to run the algorithm 924 
times for the single combination C2I (12 pulses for the 
second qubit). Hence computationally it is prohibitively 
expensive if hundreds of pulses need to be optimized. 
Second, as the pulse number increases, the optimized re- 
sult becomes more sensitive to the initial guess. Math- 
ematically it is also a challenging question to find the 
global minimum in a high-dimensional parameter space. 
For 71 + 771 — 24, we have considered two initial guesses, 
namely, an equally spaced DD sequence and a nested- 
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0.2 0.4 0.6 O.E 
optimized (4,0) 



0.2 0.4 0.6 O.E 
optimized (8,0) 



0.2 0.4 0.6 O.E 
optimized (12,0) 



0.2 0.4 0.6 O.E 
optimized (4,2) 



0.2 0.4 0.6 O.E 
optimized (8,2) 



0.2 0.4 0.6 O.E 
optimized (12,2) 



TABLE VI: Comparison between optimized two-qubit DD 
with n + m = 24 and nested-UDD(4). The meaning of the 
table columns are the same as in Table II. Noise spectrum is 
similar to that considered in Tables II and III. 



Si = L0Q{1 — U}) 


, S2 = ojQ{l - tj), S3 = 


0.5a;B(0.5 — lj) 


Pulse combination 


Pulse location 


Performance 


nested-UDD(4) 


5,10 


5.21 X 10"^ 


C24 


2,9 


2.81 X IQ-i" 


C24 


1,3,5,11 


3.31 X IQ-" 




2,4,7,8,10,12 


2.34 X IQ-" 


Si = a;e(5 - 


oj), S2 = a;0(5 — ijj),S3 


= Lje{3-uj) 


Pulse combination 


Pulse location 


Performance 


nested-UDD(4) 


5,10 


3.31 X 10"^ 


C24 


3,9 


1.42 X 10~^ 


(^24 


1,3,5,10 


1.51 X lO"'' 


(^24 


2,4,6,9,11,12 


1.35 X lO"'' 



FIG. 2: Optimized pulse timings for Si = 106(10 - cj),S2 = 
0.10(0.1 - Lu), S3 = 0.056(0.05 - Lj). The left panels are the 
results if no pulse is applied to the second qubit. The right 
panels depict the optimized DD sequences if two pulses are 
applied to the second-qubit pulses (m = 2). The grey (red) 
lines stand for timings of the pulses applied to the second 
qubit. From top to bottom, n + m = 4, 8, and 12. Note that 
in the upper right panel, the two pulses applied to the second 
qubit coincide and hence cancel each other. 

UDD(4) sequence. The associated results are shown in 
Table I VII It is seen that for the same noise spectrum, 
the performance is only slightly better than that for 
n + m = 15. It is hence possible that our result as a 
locally optimized result is still far away from the globally 
optimized result. More effective optimization algorithms 
together with more initial guesses may be the ultimate 
solution. 

It should be also noted that the total pulse number is 
physically limited. Realistic pulses are imperfect and an 
ideal instantaneous tt rotation is never possible [33| . As 
a result, the pulse-to-pulse errors may accumulate with 
increase of the pulse number (s^ . [35| and hence more 
pulses do not necessarily lead to better performance. 
Furthermore, constraint on the minimal pulse interval 
might place another intrinsic limit on the pulse number 
we could consider [sl]. For fixed total evolution time t, 
such a constraint imposes an upper bound for the pulse 
number. For a varying t, the optimized performance in 
practice will be connected with the available pulsing rate 
as well as the spectral bandwidth [s^, [l!] . 



V. CONCLUSIONS 

Using a pure dephasing model incorporating both lo- 
cal and nonlocal Gaussian noise, we have studied how the 



DD protection of two-qubit states can be better achieved 
by optimizing the pulse locations as well as the parti- 
tion of the pulse numbers allocated to each qubit. Com- 
pared with nested-UDD as a general-purpose DD scheme 
for two-qubit systems, our optimization procedure may 
improve the performance by many orders of magnitude, 
using only a few tens of instantaneous tt pulses. This 
makes it possible to use much less pulses to obtain a 
given fidelity of decoherence control in two-qubit systems. 
The price of this performance gain is the required knowl- 
edge of the noise spectrum. In addition, it is also seen 
that if the noise spectrum decays to zero slowly (e.g., 
a Lorcntzian shape), then even our optimized DD may 
not perform very well and it is only slightly better than 
nested-UDD. 

The results here also help to gain more insights into 
the issue of DD protection in two-qubit or multi-qubit 
systems. The importance of fighting against the nonlocal 
noise in two-qubit decoherence control (or more generally, 
the importance of taking into account the impact of the 
dynamics of one qubit on the environment of the other 
qubit) becomes clearer. Without the nonlocal noise, 
our optimized DD sequence simply degenerates into two 
single-qubit optimized DD sequences. In the presence of 
nonlocal noise, the optimized two-qubit DD sequence is 
quite different from two optimized single-qubit DD se- 
quences, because the pulse numbers and pulse locations 
for each qubit should be adjusted (in a somewhat sub- 
tle manner) based on both nonlocal and local noise. If 
the local noise spectrum of the two qubits is highly un- 
balanced, then most pulses should be applied to one of 
the two qubits to defeat its rapid dephasing. It is these 
specific strategies, which are not exploited by the rather 
universal nested-UDD scheme, that makes optimization 
possible. Note also that in nested-UDD schemes, the 
number of control pulses or the control order may be 
different for different control layers. But as seen from 
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one optimization result in Sec. IV-C, the performance 
of such type of asymmetric-nested-UDD is still far from 

optimized. 

There have been wide experimental interests in both 
fundamental aspects of decoherence and the DD ap- 
proach to decoherence suppression. Two types of ex- 
periments may be motivated by this work. First, for 
two-qubit systems with known noise spectrum (local and 
nonlocal), it is of immediate interest to apply optimized 
DD sequences to extend the lifetime of entangled states, 
hopefully with better efficiency. Second, one may exper- 
imentally study the nonlocal noise of a two-qubit system 
by use of single-qubit DD. That is, if two single-qubit DD 
sequences are applied on top of each other, then (assum- 
ing that the local noise of the two qubits are independent) 
only the nonlocal noise is not suppressed and hence its 
impact on two-qubit dephasing may be directly observed. 



Such type of experiments are of importance to the design 

of optimized two-qubit DD sequences and to the under- 
standing of decoherence mechanism in two-qubit systems 
embedded in a solid-state environment. 
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